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We obtain limits on the quartic neutral gauge bosons couplings Z777 and ZZ77 using LEP 2 
data published by the L3 Collaboration on the reactions e + e~ — > 777, Z77. We also obtain 95 % 



^S^ ' C. L. limits on these couplings at the future linear colliders energies. The LEP 2 data induce limits 

Oh. 
fL ' of order 10 5 for the Z777 couplings and of order 10 2 for the ZZ77 couplings, which are still 

<D ■ 

above the respective Standard Model predictions. Future e + e linear colliders may improve these 



limits by one or two orders of magnitude. 
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I. INTRODUCTION 

Neutral gauge bosons self couplings provide a window to study physics beyond the Stan- 
dard Model (SM) pj-yj]. While trilinear neutral gauge boson couplings (TNGC) ViVjVk, 
with Vi = Z,j, test the gauge structure of the SM [3J], it has been argued that quartic neu- 
tral gauge boson couplings (QNGC) VtVjVkVi may provide a connection to the mechanism of 
electroweak symmetry breaking |l( • Since the longitudinal components of the Z gauge boson 
are Goldstone bosons associated to the electroweak symmetry breaking mechanism, these 
QNGC could represent then a connection with the scalar sector of the gauge theory that has 



Q. 



become popular after the recent evidence of a new boson with a mass around 125 GeV 
However, it has been found recently in a detailed calculation of the one-loop induced decay 
mode Z — > 777, in both the SM and the 331 model, that the respective scalar contributions 
are suppressed with respect to the dominant virtual fermionic contributions [5j. This is also 
the case in the one-loop contributions to TNGC [3|, |6|]. The QNGC are induced by effective 
operators of dimension greater or equal to six and, in the SM, the QNGC are highly supp- 
resed, with the only exception of the ZZZZ vertex, because they arise at the one-loop level 
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Any deviation from the SM predictions for the QNGC will be associated to a signal 
of new physics effects [1] . 

While considerable effort has been devoted to study the TNGC, the QNGC are only 
starting to receive some attention. TNGC have been measured with an accuracy of the few 



percent level at LEP 2 [8[ and the Tevatron |9|], while QNGC are only loosely constrained at 
^EP 2 [8|. In fact, the Z777 couplings have not been bounded yet by direct measurements 
8] . In the present paper, we are interested in obtaining limits on the quartic vertices Z777 

and ZZ'-f-f coming from the LEP 2 data on the reactions e + e~ — > 777, Z77 that were used to 



get limits on the anomalous H Z7 coupling but not on the QNGC lOLIll]. We will obtain also 



95% C. L. limits on these quartic couplings at t 



re future International Linear Collider (ILC) 



and the Compact Linear Collider (CLIC) [12|, Il3| . Since there is not a published account, 



as far as we know, of the calculation of the Z777 vertices in the SM, in the Appendix we 
present a brief analysis on the connection of the Z777 form factors to the analytical results 
obtained in Ref. |5J for the branching ratio of the decay mode Z777 in both the SM and 
the 331 model. However, a similar calculation for the ZZ77 form factors in the SM is not 
available in the published literature. 



Constraints on the anomalous quartic gauge couplings ZZ77 have been studied in 77 
and Z, fusion processes at the LHC Q, in ZZ;, Z 77 production processes at future e+J 

linear colliders |6[ and from the non observation of the rare decay Z — > z/z/77 at LEP 1 [7|. 
However, constraints on the anomalous Z777 vertex are more difficult to get. In the present 
paper we find that the negative search for the reactions e + e~ — > 777, Z77 at LEP 2 by the 
L3 Collaboration may be translated into limits of order 10 -5 on the Z777 couplings and of 
order 10~ 2 on the ZZ77 couplings. We also find that sensitivity studies on these couplings 
at future e + e~ colliders may improve these limits by one or two orders of magnitude. 

The paper is organized as follows. In Section II we present the calculation of the respective 
cross sections for the processes e + e~ — > 777, Z77 and in Section III we include our results 
and conclusions. In the Appendix we give details on the connection among our quartic 
couplings Gi : 2 and the results obtained in Ref. [5| for the branching ratio of the decay mode 
Z — > 777 in the SM and the 331 model. 



II. CROSS SECTIONS 



We will use the following parameterizations for the QNGC 
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where F^ v = d^A v — d u A^ and Z^ v = d^Z y — dyZ^ are the respective gauge tensor fields for 
the photon and the Z boson. A represents the energy scale at which new physics interactions 
may appear. The respective Feynman rules for these effective vertices are thus given by, 



A 4 



A 4 



(Pi ■ P2){P2 ■ Pz)9 ap 9^ - (Pi ■ P3)PluP2p9 a p ~ (Pi ■ Pz)PluP2a9p^ + PluPlpP2pP3c 



- (pi ■ p 2 )(Pl ■ Pz)g a »9p V + (P2 • V?)P\aP\v9pp. ~ (P2 • P3)PluPl P g ap + (>2 • Pz)P\vP2<x9pp 



+2(p 2 ■ p 3 )plpP2pgau ~ (Pi ■ P3)PlaP2 P gpu ~ PlaPluP2pP3p 



(3) 



and 



ic 



Sc^A 2 



Aa g c 



(pi-P2)g*"-rfj% 



+ a c 



-i Q -n^ _l_ ^ r , a \nt iu 



(pip p 2 + m)gr + (pi ■ p 2 )(<r<r + <r <? 



lia vfi 1 va iif}\ 



- pftp&T + VW ? ) - didsT + Pi9^) ] } • 



(4) 



All the couplings G\p and a 0>c are CP conserving and within the SM all of them vanish at 
tree level. As far as we know, the ao, c have not been computed explicitly in the SM, whereas 
the couplings G\ t 2 can be extracted directly from the recent calculation performed in the 
SM and the 331 model [16J for the branching ratio of the rare decay mode Z —¥ 777 5J. 
Since these authors did not use explicitly the parametrization given in Eqs. (JTJ) and (121), we 
have included in the Appendix the connection of the G\^ couplings to the results obtained 
for the Z — y 777 decay in Ref. [5]. These form factors are dominated by the fermionic 
virtual contributions and they are essentially the same in both the SM and the 331 model, 
but unfortunately with rather low values, 1.63 x 10~ 10 and 1.33 x 10~ 10 , respectively. 

In Figures [1] and [2] we present the contributions of the effective interactions given in Eqs. 
© and (HD to the processes e + e~ — > 777 and e + e~ — > Z77. The SM contributions to these 

nn 

processes occur via t-channel diagrams involving initial-state radiation [2, 6]. The respective 
SM cross sections have been computed by Stirling and Werthenbach for CM energies greater 
than 200 GeV [6j. They are of order of few femtobarns. According to this result, in 
order to reduce the contributions due to initial-state radiation in these reactions, the L3 
Collaboration introduced cuts on the photon energies and their polar angles, E 1 > 5 GeV 
and I cos^l < 0.97 [10]. Events from e + e~ — > 777, Z77 processes were selected by requiring 
the photon candidates to lay in the central region of the detector with | cos# 7 | < 0.8 and a 
total CM electromagnetic energy large than a/s/2. In this case, the L3 Collaboration was 
interested in getting limits on the anomalous Higgs couplings H Z7 and if 77. However, using 
their data we are able to get also limits on the G\ t 2 and ao, c couplings: G1/A 4 < 1.2 x 10~ 5 , 
G 2 /A 4 < 9.4 x 10" 6 , a /A 2 < 5.9 x 10" 3 and a c /A 2 < 1.6 x 10~ 2 . The latter limits should be 
compared with the more stringent bounds obtained by the L3 Collaboration from a direct 
search of Z77 events at LEP 2 energies: [-0.08, 0.021] and [-0.029, 0.039], respectively [sj. 

The expressions for the respective cross sections induced by the effective vertices given 
in Eqs. ([3]) and (j4j) are given by 
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In Eq. ([()]) , the first term comes from the Feynman diagrams shown in Fig. [2] for the 
exchanged photon and the second one comes from the exchanged Z boson. We did not 
include the contribution coming from the interference of the two diagrams because we will 
get limits on the form factors one at the time. 

III. RESULTS AND CONCLUSIONS 



Using 



the numerical values sin 2 9w 



Xw 



0.2314, M z = 91.18 GeV and T z = 2.49 



GeV [17] . we obtain the cross sections for the processes e + e~ — > 777, Z77 as functions of 
the CM energy and the Gi,2, a 0jC couplings. We have also implemented in our calculation 
the cut used by the L3 Collaboration on the CM energy and the photon polar angle in order 
to suppress the SM contributions associated to initial-state radiation. In Fig. [3] we depict 
the sensitivity limits at 95 % C.L. for the G\^ couplings for CM energies of 500 GeV and 
1000 GeV and we have taken the G\$, couplings one at the time. The respective combined 
limits contours are shown in Fig. HJ On the other hand, in order to get sensitivity limits 
and the respective limits contours for the ao, c couplings we have set to zero the contribution 
associated to the G\$ couplings in Eq. fl6]). The respective limits are given in Figures [5] and 
Oalso for two different values of the energy of the ILC and CLIC acelerators. 

In conclusion, we have obtained limits on the quartic couplings Z777 and ZZ77 at LEP 
2 energies by using published L3 data for the reactions e + e~ — > 777, Z77. However, our 
limits for the Z777 couplings are five orders of magnitude above the expected predictions 
for the SM included in our Appendix. In the case of the ZZ77 couplings, our limits obtained 
from the LEP 2 data on the reaction e + e~ — > Z77 are still one order of magnitude above the 
best limits obtained in LEP collider [8]. On the other hand, these limits may be improved 
by one or two orders of magnitude at future linear colliders (Figure EH). Similar limits have 
been obtained for the a 0jC couplings from the process e + e~ — > z/z/77 [18] and through effects 



induced by the polarization of the Z gauge boson and initial state radiation in the process 
e+e^ -> Z77 [12]. 
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Appendix A: Z — > 777 at one-loop level 
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The branching ratio for the decay Z — > 777 was obtained in Ref. [5| for the SM and the 
331 model to one- loop order. In this appendix we present some details of the calculation 
that leads to the values obtained for the quartic coupling constant Gj/A 4 , % = 1,2, defined 
in Eq. ([1]). In Ref. J5( this calculation was performed in terms of the form factors F Zi 
which identify the fermionic, vectorial, and scalar components of the respective Feynman 
amplitudes, 



F^ + F'+Fl, (Al) 



'Zi = r z% -r r Zi 



with 



F l z{ 2 = zZ 9i/2fl = Fl+Fg+F Zt . (A2) 

F=q,Q,l 

1/2 

The fermionic contribution F z [ is expresed in terms of the light quarks q = u,d } s, c, b, t, 
the new heavy quarks of the 331 model Q = D, S, T, and the leptones I = e, fi, r, with 

gf /2 = -N F Qlg F z , (A3) 

where Q Fi is the respective fermionic charge and g z correspond to the fermionic coupling 
constants to the neutral gauge boson Z and Z'. The vectorial form factor F z _ is splited in 
the contributions of the charged W boson and the respective degenerated bileptons Y + and 
Y ++ , and their correspondent pseudoscalars 



Fk = Ff l+ Fl 

= [2cl.fl - c 2W f°y] + [(7 - 3As 2 w )(fl - f z J)\ . (A4) 

Finally, the scalar form factor F z _ is given in terms of the charged Higgs bosons H = hf , 
h + h + hi d ++ d ++ d ++ 



r Zi — r Zi 



2(9 - h2s 2 w )f h Zi . (A5) 



The form factors f%. are given in Ref. 19j, and f z .' ' ' 'in the appendix of Ref. [5(. 
The relation between the effective form factors G,/A 4 and the Fz form factors calculated 

n 

in Ref. [5J is given by the integral expression 
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where the x and y are the kinematical variables asociated to the Z — y 777 decay amplitude 
J5|], and p = e/sw = V^ira/sw- 

In the Tables Hl lVIII we show the results obtained for each one of the contributions of the 
Gj/A 4 form factors for the SM and the 331 model. 

The branching ratio predicted for the Z — y 777 decay in these models is obtained with 
the relation 



Bi(Z —y 777) 



T(Z ->■ 777) 



and the decay width [15|, [2 



r(Z — y 777) 



Ml (2G\ + 3G 2 2 - 3GiG 2 ) 



^2960vr 3 
and the experimental total Z width 17|] is 



A 8 



(A7) 



(A8) 



2.4952 GeV. 



(A9) 



In Table IVIIII we give the predictions obtained for the G;/A 4 [GeV 4 ] form factor from 
the branching ratio (A7). 



TABLE I: Predictions of the 331 model to the Gi^/A 4 [GeV 4 ] couplings and the respective 
BR(Z — > 777) from the fermionic, bosonic and scalar sectors. BR with Eq. (|A8|) , 



Sector |d/A 4 | |G 2 /A 4 | BR 

Fermions 6.06xlCT 10 8.66xlCT 10 1.44xlCT 8 

Gauge Bosons 1.41xl(r 12 1.75xl0 _n 8.65xl(T 12 

Scalar 6.24xl0~ 14 1.27xl0" 13 3.30xl(T 16 

Fermions-Gauge Bosons 2.65X1CT 11 1.12xlCT 10 3.07xlCT 10 

Fermions-Scalar 6.10xl(T 12 1.15X10 -11 2.66xl(T 12 

Gauge Bosons-Scalar 1.82xl(r 13 4.37xl(T 13 4.08xl(r 15 

Total 6.06xl(T 10 8.74xl(r 10 1.46xl(T 8 



TABLE II: Predictions of the 331 model to the Gi^/A 4 [GeV 4 ] couplings and the respective 
BR(Z — > 777) from fermionic subsectors. BR with Eq. (|A8jl . 



Fermions |Gi/A 4 | |G 2 /A 4 | BR 

Quarks 4.75xl(T 10 6.73xl(T 10 8.67xl(T 9 

Leptones 1.42xl(r 10 2.03xl(T 10 7.89xl(T 10 

Quarks-Leptons 3.50xl(r 10 5.06xl(r 10 4.91xl(T 9 
Total 6.06xl(T 10 8.66xl(T 10 1.44xl(T 8 



TABLE III: Predictions of the 331 model to the Gi^/A 4 [GeV 4 ] couplings and the respective 
BR(Z — > 777) from quarks subsectors. BR with Eq. ()A8|h 



Quarks |Gi/A 4 | |G 2 /A 4 | BR 

SM quarks 4.73xl0" 10 6.70xl(T 10 8.59xl0" 9 

Exotic quarks 2.35xl(T 12 5.15xl(T 12 5.53xl(T 13 
SM-exotic quarks 3.35xKT n 6.28xKT n 7.91xl(r n 
Total 4.75xl(T 10 6.73xl(T 10 8.67xl(T 9 



TABLE IV: Predictions of the 331 model to the Gi^/A 4 [GeV 4 ] couplings and the respective 
BR(Z — > 777) from gauge bosons. BR with Eq. (1A8D . 



Bosons |Gi/A 4 | |G 2 /A 4 | BR 

W 1.42xl(T 12 1.76xl0 _11 8.75xl(T 12 

Bileptons 4.44xl(T 14 1.28xl(T 13 3.67xl(T 16 
W-Bileptons 1.66xl(T 13 1.79xl(T 12 8.94xl(T 14 
Total 1.41 xlO -12 1.75XKT 11 8.65xl(T 12 
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TABLE V: Predictions of the SM to the Gi, 2 /A 4 [GeV~ 4 ] couplings and the respective BR(Z 
777). BR with Eq. (|A"8J) . 



SM |Gi/A 4 | |G 2 /A 4 | BR 

Quarks 4.73xl(T 10 6.70xl0" 10 8.59xl0" 9 

Leptons 1.42xHT 10 2.03xl(r 10 7.89xl(T 10 

W 1A2x10- 12 1.76xl0 -11 8.75XKT 12 
Quarks-Leptons 3.49xl(T 10 5.05xl(T 10 4.89xl(T 9 

Quarks-W 2.31xl0 _n l.OOxlO -10 2.46xl(T 10 

Leptons-W 1.19X10 -11 5.16X10 -11 6.55xHT n 

Total 6.04xl(T 10 8.71 xlO -10 1.45xl0 -8 



TABLE VI: Predictions of the SM to the Gi, 2 /A 4 [GeF~ 4 ] couplings and the respective BR(Z 
777) from quarks. BR with Eq. (|A8[) . 



SM quarks |d/A 4 | |G 2 /A 4 | BR 

3.11xl(T 10 4.43xl(T 10 3.76xl(T 9 
2.58xl(T 10 3.18xl(T 10 1.94xl0 -9 

t 1.85xHT 14 2.11xl(r 14 8.65xl(T 18 



u 



c 



7.02xl0 _n l.OOxlO -10 1.92xl0 -10 
7.50xl0 _11 l.OTxlO -10 2.19xl0 -10 
1.86xlO _n 2.43X10 -11 1.13xl0 -11 
Interference 2.23 xlO -10 3.60xl(T 10 2.52xl0 -9 
Total 4.73xHT 10 6.70xl(r 10 8.59xl0 -9 



11 



TABLE VII: Predictions of the SM to the Gi^/A 4 [GeV~ 4 ] couplings and the respective BR(Z 
777) from leptons. BR with Eq. (|A8|) . 



Leptons |Gi/A 4 | |G 2 /A 4 | BR 

e 6.88X10" 11 9.80X10" 11 1.84xl0 -10 

fj. 7.39xl0 -11 1.06xHT 10 2.15xl(T 10 

r 4.05xl0~ n 5.10xl0 -11 4.98XHT 11 

Interf. 9.12xl0 -11 1.33xl(T 10 3.39xl(T 10 
Total 1.42xl0 -10 2.03xl(T 10 7.89xl(T 10 



TABLE VIII: Values of G ii2 /A 4 [GeV~ 4 ] as function of BR according to Eq. (16) of Ref. jlj]. We 
include the PDG 2012 limit for BR(Z ->• 777) [lj. 



BR |Gi/A 4 | |G 2 /A 4 | 

1(T 5 (PDG) 2.22xl(T 8 1.81xl(T 8 

5.41xl(T 10 (SM) 1.63xl(T 10 1.33X10 -10 

5.26xl(T 10 (3-3-1) 1.61xl(T 10 1.31xl0 -10 
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FIG. 1: Feynman diagram for the process e + e — > 777 induced by the effective vertex Z777. 
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FIG. 2: Feynman diagrams for the process e + e —> Zjj induced by the effective vertices ZZ'yy 
and Z777. 
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FIG. 3: Sensitivity limits at 95 % C.L. for the couplings -^-[GeV^ 4 ] as function of the integrated 
luminosity for two ILC/CLIC CM energies. 
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FIG. 4: Contours limits at 95 % C. L. in the G1-G2 plane for the process e + e — > 777 for two 
values of the cross section, 40 and 100 fb and yfs = 500 GeV . 
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